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1 Introduction 

Let Xi, . . . , Xn be i.i.d. observations, where Xi = Yi + Zi and Yi and Zi are independent random 
variables. Assume that the unobservable Y^s have distribution function F and density /, and 
that the Z's have a known density k. Note that g equals k * f, where * denotes convolution. 
The deconvolution problem is the problem of estimating the density / from the observations 
Xi from the convolution density g. 

A well known estimator of f{x) is based on Fourier inversion and kernel smoothing. Let 
w denote a kernel function and h > a. bandwidth. The kernel type estimator fnh{x) of the 
density / at the point x is defined as 
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with 

„,.(,.) ^ f e-'",s. 

27r y_oo <Pk{s/h) 

Here (pemp denotes the empirical characteristic function of the sample, i.e. 4>emp(t) = ;^ J2]j=i ' 
and (pw and (pk denote the characteristic functions of w and k respectively. Note that, even 
though has the form of an ordinary kernel density estimator, because of the dependence 
of Vh on the bandwidth h it is different. Kernel type estimators for the density / and its dis- 
tribution function F have been studied by many authors. Relatively recent papers are Zhang 
(1990), Fan (1991a,b), Fan and Liu (1997), Van Es and Kok (1998), Cator (2001), Van Es and 
Uh (2001), and Delaigle and Gijbels (2002). For an introduction see Wand and Jones (1995). 
This paper covers a chapter in Uh (2003). 

The expectation of the estimator (jl.lj) has a familiar form. We have, see for instance 
Stefanski and Carroll (1990), 

EUix) = E^w[^^). (1.2) 

Indeed, this expectation is equal to the expectation of an ordinary kernel density estimator of / 
based on observations Yj from /. Expansions of (|1.2p for h ^ are standard in kernel density 
estimation theory and are hence omitted here. See for instance Wand and Jones (1995). 

Deconvolution problems are usually divided in two groups, ordinary smooth deconvolution 
problems, where the rate of decay to zero at infinity and minus infinity of the characteris- 
tic function 0^ is algebraic, and super smooth deconvolution problems, where it is essentially 
exponential. This rate of decay, and hence the smoothness of the known density /c, has a 
tremendous influence on the variance of the estimator, see for instance Fan (1991) or Cator 
(2001). By (jl.2|) it is clear that the expectation is not affected. The general picture is that 
with increasing smoothness of k the estimation problem becomes harder and the the optimal 
rates become slower. 

Our aim is to derive classical central limit type theorems for these kernel type deconvolution 
estimators. For ordinary smooth deconvolution this has first been achieved in Fan (1991) and 
extended in Fan and Liu (1997). The limit behaviour in this case is essentially equal to that 
of a kernel estimator of a higher order derivative of a density. In some specific deconvolution 
problems this is evident from relatively simple inversion formulas, cf. Van Es and Kok (1998). 
For instance, for generalized gamma deconvolution where k is the density of Aii^i + A2-E2 + ■ ■ ■ + 
XmEm, with Ai > 1, . . . , Am > and Ei, . . . , independent standard exponential random 
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variables we have 

V^J^ifnhix) - E U{x)) ^ N{0, a'), (1.3) 

where = (Ai . . . A^)^ / w^"^\vydvg{x). This result is typical for ordinary smooth deconvo- 
lution, in the sense of a rate of convergence that is algebraic in h. 

Asymptotic normality of fnh{x) in super smooth deconvolution problems has been derived 
by Zhang (1991), Fan (1991b) and Fan and Liu (1997). Under suitable conditions their theorems 
state 

^flllM^lIluMZNiO,!), (1.4) 

where Znj = j['"h(^^-^^, j = l,---n and either = ^Yl]=i^n.j equals the sample 

variance of So the estimator is studentized in some respect. The asymptotic 

variance is not clear. Van Es and Uh (2001) derive a central limit type theorem like ()1.3p for 
super smooth deconvolution, where the asymptotic variance is clear and the normalisation is 
deterministic. Their result is given in Theorem 11.11 below. 

Condition W 

Let 0^ be real valued, symmetric and have support [—1, 1]. Let 0tu(O) = 1, and let 
0^(1 — t) = At" + o(t") as it J, 0, for some constants A and a > 0. 

Condition K 

Assume that (pk has exponentially decreasing tails, i.e. 0fc(t) ~ C|t|'*'°e^'*'^'''^, as |t| oo, 
for some A > 1,/i > 0, Aq, and some real constant C. Furthermore assume (pkit) ^ for all t. 

Note that Condition K excludes the Cauchy distribution and all other distributions for which 
the tail of the characteristic function decreases more slowly than e~'*'. 

Theorem 1.1. Assume Condition W, Condition K and EX^ < oo. Then, as n oo and 

^ iUix) - E U{x)) ^ N{0, il(^/A)2+2«^r^^ ^ i))2)_ ^15) 

^A(i+a)+Ao-ig;:^ ivr^ 

Surprisingly, the asymptotic variance is distribution free, in the sense that it does not depend 
on / or X. The condition A > 1 was needed to ensure that remainder terms in the proof of this 
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theorem are asymptotically negligible. Note also the condition that the second moment of the 
observations is finite. 

By studying deconvolution problems where the known distribution is a symmetric stable 
distribution we will investigate the asymptotic behavior of the kernel deconvolution estimators 
if the conditions of Theorem 11.11 are not satisfied. 

2 Deconvolution for symmetric stable densities 

Consider deconvolution for symmetric stable densities k which have characteristic function 

(l),^{t) = e-l*l^/^, /X > and < A < 2. (2.1) 

The condition < A < 2 is necessary to ensure that is a density, cf. Chung (1974). Hence 
the normal distribution is, in some sense, extreme. Note that for A equal to one is a Cauchy 
density. The only symmetric stable distribution with finite second moment is the normal 
distribution for which A equals two. This implies that, unless A equals two, the second moment 
of the observations will be infinite. Hence the normal distribution is the only symmetric stable 
distribution for which Theorem 11.11 applies. We will derive a limit behavior, similar to that 
described by Theorem ll.il for cases where A is larger than one. Of even more interest are the 
cases where A is equal to one, i.e. Cauchy deconvolution, or smaller than one. It turns out 
that, while crossing the Cauchy boundary, a different limit behavior appears. 

For simplicity we only consider the sine kernel, defined by 

w{x) = sinx/(7rx) {4>w(t) = /[_i,i](i(:)). (2.2) 

Results for a more general class of kernels w are given in Uh (2003). 

First we give a heuristic derivation of the results which are rigorously proved in Section IHl 
Note that the estimator fnh can be rewritten as 

/-(^•)^^^i:fcos(.(izf:)),w«V,.,,. p.3) 

Let S denote a random variable, independent of the Xj, having probability density fs given by 

^""^^ c{h) 
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on the interval [0, 1], where the normahzation constant c{h) is given 
Then (j2.H|l is equal to 



fnhix) = ^ Ve(cos((S//i)(X, -x))|X,) = ^ VE5Cos((S//i)(X, -x)). (2.4) 



It turns out that the asymptotics are greatly determined by the asymptotics of the distribu- 
tion of the random variable S. Let us first consider its expectation and variance. The following 
lemma gives expansions of the normalization constant, the expectation of S and the variance 
of S. Its proof is given in Section HI 

Lemma 2.1. For < A < 2 and h ^ we have 

c{h) = (^^/^^ + 0(/^2^))e(l/'^)'/^ (2.5) 

= l-^h^ + o{h^), (2.6) 
A 

Var5 = ^h^^ + o(h^^). (2.7) 
These expansions suggest to normalize S as follows. Write 

E„ ^ (2.8, 
The density function, say, of En is given by 

jEr^yv) - ^ ^^^^ e i[-(A^)/^-^o]l^^J• 

By Taylor expansion and Lemma lTTl it converges uniformly on bounded intervals to e'"I(_oo,o]('^)- 
This implies that En converges in distribution to —E where E denotes a standard exponential 
random variable. 

For the terms in ()2.4|) we have 

Escos{{S/h){Xj-x)) = Escos{{{l + S -l)/h){Xj-x)) 
= E s cos{{l/h){Xj - x)) cos{{{S - l)/h){Xj - x)) 

-Essm{{l/h){Xj-x))sm{{{S ~l)/h){Xj-x)) 
= cos (^^1^)Es„ cos h'-'En{X, - x)) (2.9) 

- sin (^^)Ee„ sin h^-'En{X, - x)) . (2.10) 
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It now becomes apparent that we may expect different asymptotics in the cases < A < 1, 
A = 1 and 1 < A < 2. In these cases the factor h^^^ in (j2.9p and (|2.1(jp diverges to infinity, 
equals one and vanishes. 

The next three theorems establish asymptotic normality for 1 < A < 2, i.e. for the symmetric 
stable densities whose characteristic function decreases more rapidly than the characteristic 
function of the Cauchy distribution, for Cauchy deconvolution, and for 1/3 < A < 1, i.e. for 
the symmetric stable densities whose characteristic function decreases more slowly than the 
characteristic function of the Cauchy distribution. 



Theorem 2.1. Let w be the sine kernel \2.2^) . // 1 < A < 2 then, as n ^ oo and h 0, we 
have 

— {Uh{x) - E^,(x)) ^ iV(o, —A^l|\f). (2.11) 



fi^-if.'^ 27r^ 



Theorem 2.2. Letw he the sine kernel If X equals one, i.e. Cauchy deconvolution, then, 

as n ^ CO and h 0, we have 

V^e-^Afnh{x) - EU{x)) ^ N{0,a^), (2.12) 

with 

,2 _ 1 /■ 



Theorem 2.3. Let w be the sine kernel h2.i^) . If 1/3 < A < 1 then, as n ^ 00, h — » and 

nh — s> 00, we have 

(/„;,(x)-E/„,(x))^iV(0,a2), (2.14) 



/^(A-l)/2g(l//i)VM 

with 



a' = ^^gix). (2.15) 



The global picture we see from these three theorems is that for 1/3 < A < 1, apart from the 
exponential rate of convergence, the asymptotic variance resembles the asymptotic variance of 
a kernel density estimator, in the sense that it depends on the value of g at the point x, as 
in ()1.3|) . This is typical for smooth deconvolution problems, though the rate of the variance 
is exponential in h and not algebraic. For Cauchy deconvolution we see that the asymptotic 
variance depends globally on g. For 1 < A < 2, the estimator is asymptotically distribution free. 
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It shares the asymptotics of Theorem ll.il even though the second moment of the observations is 
infinite for 1 < A < 2. Concluding we see that the restriction A > 1 in Theorem II .11 is essential 
and that the finite second moment condition might not be. Crossing the Cauchy boundary we 
get different asymptotics. 

3 Proofs 

3.1 Basic lemma 

Lemma 3.1. Let = {X — x)/h mod 27r. As n ^ oo and h ^ 0, 

{X,Y,) ^ {X,U), (3.1) 

where U is uniformly distributed on the interval [0,27r]. Moreover X and U are independent. 

Assume < A < 1. Let z he a hounded periodic function with period 27r and let w he a 
continuous and integrahle function such that w is monotone in the tails. Then, as n ^ oo and 
h ^ 0, we have 

h^-^E (^z(J^^^^'^w{h^-\X, - x))'^ ^ ^ j\{u)du j w{u)dug{x). (3.2) 

Proof 

Note that the density g = k* f of X is continuous and bounded. For — oo <m<oo, 0<?/<27r 
and M < u, we have by a Riemann sum approximation of the integral of g over the interval 

[M,u], 

f{2ni+y)h+x 

P{M<X<u,Y,,<y)= Yl / . ' 9{t)dt + 0{h) 

i:M<(2m+y)h+x<u >^27rj/i+x 

E yhg{^^,h) + 0{h) = ^ 27Thg{^,,,) + 0{h) 

i:M<(2ni+y)h+x<u i:M<{2m+y)h+x<u 

= ^ r g{u) du + 0(1) = ^ {G{u) - G{M)) + o(l), 

where E,i,h is a point on the interval [2z7r/i + x, 2z7r/i + yh + x\ C [2z7r/i + x, 2{i + l)Tih + x]. 

For arbitrary e > choose M(e) such that G{M{e)) = P{X < M(e)) < | e and no{e) such 
that for n > no{e) we have |P(M(e) < X < u,Yh < y) - g {G{u) - G(M(e)))| < | e. Then, 
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for n > 'n.o(e), 

\P{X <u,Yn<y)-^G{u)\ 

< \P{X <u,Yn<y)- P{M{e) < X < u,Yh < y)\ 

+ |P(M(e) <X<u,Y,<y)~^ - G(M(e)))| 

+ \^{G{u)-G{M{e)))-^G{u)\ 

< P{X < M(e), < y) + i e + G{M{e)) < e. 

Hence limP(X < u,Yh < y) = ^ G{u), which proves (j3.ip . 
To prove ()3.2|) write 

h^-^Ez^^^^^wih'^-'iX.-x)) 



I z{v)w{h^v)g{x + hv)dv 

f27r 



z{t + 2ni)w{h^{t + 2ni))g{x + h{t + 2TTi))dt 



— I z{t)p{t)dt, 







where 



pit) = 2nh^ w{h^{t + 2m))g{x + h{t + 2Tci)). (3.3) 
With ^i(t) = h^{t + 27ii) for i e Z, we have 

Pit) = 27r/i" Yl ^(Ut))9{^ + h'-'m) 
= gix) 2nh^ Yl + 2^^^ Y.'^iUt)){9{x + h'~%it)) - gix)). (3.4) 



Let M > be such that \w\ is increasing on (— oo, —M] and decreasing on [M, oo). Then, 
for < t < 2tt, we have 

\wiu)\du + oil) (3.5) 

and 



2nh' Y Mm)\<^^h' Y Hm)\+oih')< 



\wiu)\du + oil). (3.6) 



i:ii(t)>M i-ii{0)>M 



The convergence of the sum to the integral follows from the approximation from below of w by 
a step function and the dominated convergence theorem. The o(l) terms in (|3.5p and (|3.6|) can 
be chosen such that they do not depend on t. 
Moreover, note that 

/M 
\w{u)\du, (3.7) 

uniformly for t in [0, 27r], by the continuity of w and Riemann sum approximation. 
By the bounds ()3.5|) and ()3.6|) . and the uniform convergence in ()3.7|) . one can show 

/oo 
\w{u)\du, (3.8) 

uniformly for t in [0, 27r]. Similarly one can show 

/oo 
w{u)du, (3.9) 

uniformly for t in [0, 27r], which implies that the first term in ()3.4|) converges to g{x) w{u)du, 
uniformly for t in [0, 27r]. Since is uniformly continuous we have g{x + /i^^'^Ct(t)) — (7(2;) — 0, 
uniformly for i and t with — M < ,^j(t) < M. Using ()3.5j) . ()3.(ij) and ()3.8|1 one can show that 
the second term in ()3.4|) vanishes, uniformly for t in [0, 27r]. Hence p{t) (^(x) f^^w{u)du, 
uniformly in t. 
Finally we get 

h^-'Ez(^^^^)w(h^-\Xi - x)) = — r zit)pit)dt 
\ h J 2tt Jq 

/oo 
w{u)du g{x), (3.10) 
-00 

which completes the proof of the Lemma. □ 
3.2 Proofs of Theorems [O, EH and 1^ 

We can derive a bound on the error in substituting —E for in the terms ()2.9p and ()2.10|) . 
The proof is given in Section EJ 

Lemma 3.2. //O < A < 2, as n 00 and h 0, we have almost surely 

E E„ cos h^'^EniXj - a;)) - E s cos ( - ^ h^~^E{X, - x))\= 0{h^) 



1 

2^ 



271 







and 



E sin /i^-^E„(X, - x)) - E sin ( - ^ h^-'E{X^ ' ^^)\ = ^^^^^ 



,A ' V ^ V A 

where E is a standard exponential random variable. 



We can now approximate fnh{x) — E fnh{x). We have 



= ^J2{^scos{{S/h){Xj-x))-EEscos{{S/h){Xj-x)) 



c{K)_ 
Tinh 



n 



cos 1 ^^)e^„ cos h'-'E^iX^ - x] 



sm 



Xi — a; 



Es„sin(^/i"-^^„(X,-x; 



- E ( cos (^^) E cos h'-'E^iX, - x] 



c{h) 



Tinh 



+ E (sin (^^)E^„sin [^h'-'E^X, - x)))) 
cos ( ^LZ^^E^cos ( - ^ - a;; 



— sm 



Xj-x 



EEsm(^-^h^-^E{Xj-x)^^ 



_ E ( cos (^^) E B cos ( - ^ h^-^E{X, - x] 
+ E ( sin (^^) E ^ sin ( - ^ h^'^E{X^ - x) 



(3.11) 

The order of the remainder term follows from the fact that it is equal to an average of inde- 
pendent terms, each of which is equal to the sum of 

'X — X'' 



c(h) ( f Xj — X 
' cos ' 



nh 
and 

c{h) 



h 



Ee,, cos h^'^En{X,-x)) -cos ^ 



h 



Eecos (^-jh^~^E{Xj-x) 



sm 



Xi — X 



Ee„ sin ( ^ h^-'En{X,-x)]-sm ' 



X~ — X 



EESin[-^h^-^ES{Xj-x) 



Tlh \ \ h J ^" \X "V ^ V \ h 

minus their expectations. The variances of these terms are of order c{hY0{h'^^)/{7c'^h'^) by 
LemmaEH which is of order Oih^^-^e^^^^^^^/^") by (Q. Hence the variance of the average is 
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of order 0(^^^'*'~^e^^^/^-'^/'^), which yields the order of the remainder term ()3.1H) by the Markov 
inequahty. 

A straightforward computation yields 

and 



where 



= -— — — — and w;2('u) - 



+ /i^M^ A^ + /i^M^ 

Note that wi and W2 are continuous bounded functions with wi{uY + W2{uY = wi{u). Note 
also that W2 is not integrable, so we can not apply Lemma (3.11 directly for w = W2- However, 
since W2 is integrable for a > 1, it turns out that we can circumvent this problem. 
Define the random variables Vnj as 

= cos (^^)w;i(/i^-i(X, - x)) - Sin (^^^i^'^w2{h^-\X, - x)) 

= cos{Yh)wiih^'\Xj - x)) - smiYh)w2{h''-\Xj - x)). (3.12) 



Then 

U{X) - E Uix) = 1 f^iVn, - E + Op /,2A-lg{l/h)VMA . (3.13) 

j=i 



To prove our three theorems we will check the Lyapounov condition for ^ X]j'=i(^",j ^ E\4j) 
to be asymptotically normal, i.e. for some 5 > we have to check 

1 2+5 

0. (3.14) 



E \ Vn,j — EKiji 



n'5/2(Var(V;j))'+^/' 

We will check this condition for 5 equal to two. Note that by the inequality |a + < 2'''(|a|''^ + 
\h\P),p > we have E(V;,i - EK,,)' =< K^V^^ + (EK,,,)')- 

For 1 < A < 2 the factor h^~^ vanishes. Hence, by Lemma 13.11 we have {h^^^X^ Yh) ^ 
(0, U). Since we are dealing with bounded continuous functions of {h^~^X, Yh) we also have 

E = E (cos(n)^i(/i^-i(X, - x)) - sm{Yh)w2{h^~\Xj - x))) 
E (cos(f/)w;i(0) - sin(f/)w;2(0)) = 

and 
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EK' = E (cos(n)w;i(/i^-i(X, - x)) - sm{Yn)w2ih^-\X, - x)))^ 

E (cos(f/)wi(0) - sm{U)w2{0)y = E cos(f/)^ = |. 

8 

The asymptotic variance is given by 

Var(V;,,) = ^Vlj - (EV;,,)2 E(cos(?7)t/7i(0) - sin(f/)w;2(0))2 = E cos{U)^ = i. 

Let us check ()H.14j) with 6 equal to two. Indeed we have 

E|K,-EV;,/ _ 0(1) . ..... 

n(Var(V;,,))2 n(i + o(l))2 " ^ ' ^ 

This shows that ^ Z]"=i(Ki,j - E Kj) and ^ ^ I]"=i(K,i-E V;,i) are asymptotically normally 
distributed. The asymptotic variance of ^ X]r=i(^n,i ~ E\4.j) is given by 

V-(^-i:(V;. -EK,)) = - ^ Var(K,i) ~ i -L^ e^(V^)V.. (3.16) 

by Lemma [2. 11 

Now consider Cauchy deconvolution where A equals one. By Lemma IXTl since we are dealing 
with bounded continuous functions of {X,Yh), we have 

EVnj = E {cos{Yh)w,{X -x)- sm{Yh)w2{X - x)) 

E {cos{U)wi{X - x) - sm{U)w2iX - x)) = 

and 

EV^j = E (cos(n>i(X - x) - sin(r;,)u;2(X - x)Y 
E {cos{U)wi{X - x) - sm{U)w2{X - x))\ 

The asymptotic variance of V^j is given by 

Var(K,,)=EK',-(EK,,)' 
^ E {cos{U)wi{X -x)- sm{U)w2{X - x))^ 
= E cos{U)^Ewi{X -xY + E sm{U)^Ew2{X - x)^ 
= - Et(;i(X - x) = - / — - — — -Tg{u)du. 



2 ' 2 J l + i2^{u-xy 

As above this shows that ()3.14|) is satisfied for 6 equal to two. Hence ^ X]j=i(^,i ~ E V^j) and 



n Sj=i(^,« ~ El^ j) are asymptotically normally distributed. The asymptotic variance of 
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^^Er=i(K,-EK,) is given by 
'c{h) 1 ^ 



^ Var(K,i) ~ - e^/(-'^) ^ / — ^ 9{u)du. (3.17) 



4-Z 



n TT^/i^ ' n 2tt'^ J 1 + /^^(u — x)^ 

by Lemma [2. 11 

Note that, if < A < 1, the factor /i'^"^ diverges to infinity. In this case we have 

EKf. = E (^)e cos (^^) Vl)^-'sin [^^y~'w,{h'-\X,-x)yw,{h'~\X,-x)) 

Since, for / = 0, 1, . . . , 4, the functions w{w2~'' are integrable and monotone in the tails, by 
Lemma [3. II we get 

= 0(/ii-^). (3.18) 

By a similar argument we have ^V^j = 0{h'^~^), and hence (El/„j)^ < {EV^j)"^ = 0(/i^"^^) = 
Oih'-^). 

Next let us consider Var(V^j) = EVj^j — (EVnjY- Using the inequality above for p = 3/2, 
the fact that w'^^'^ is integrable, and Lemma IXTl we get E|\4 .,|^/^ = 0{h^^^). By the Jensen 
inequality we have |E V;,,f < (E |Kj|)'/' < E |V;,,f /2, and so {EVnjf < ((E |Kjf /')'/')' = 
0(/i(4/3)(i-A^^_ Moreover, by (jS3) we get 

E(cos ~^ )V(/^^~'(Xj - a:))' = ^h^-^ J wi{ufdug{x) + o{h^-^). 
Similarly we have 

E(sin(^^^)'t/72(/i^"'(X, -x))2 = l/ii-^y u'2(M)'t^M^7(x) + o(/i^-^). 

and 

-)wi(/l^"^(Xj - x))w2(/i^~^(Xi - X)) 



/,A-iE COS f^V^") sin 



/I 7 \ h 
sin('u) cos(-u)(i'U / Wi(-u)w2('u)(i'U(7(x) + o(/i^^^) = o(/i^^'*') 







2 

Hence 

' /i\ \\ . ( i ^ 



EV;2. = E( COS 



(cos l^^^^j-^^Wiih^'-^Xj - x)) - sin (^lhL_^^w2{h^-\X, ~ x))y 



^ h}-^ j {wl{u) + wl{u))dug{x) + o{h^-^) = i /i^-^ / wiiu)dug{x) + o{h^-^) 
^^nh'-'g{x) + o{h'-'), 
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and 



Var(K,,,) = ^-vr/ii-^(7(x) + o(/ii-^) +0(/i(^/=')(i-^)) = \- ^h^'^ g{x) + o{h^-''). (3.19) 



Finally we check ()3.14p with 5 equal to two. Indeed we have 

E|V;,,-EK,,f ^ 0{h^-^) ^ //i^x 

nVar(K,i))2 n(|^ tt/^i-^ ^(x) + o(/ii-^))2 U/iy 



0. (3.20) 



This shows that \ Zl"=i(Kj -E Kj) and ^ ^ Zl"=i(Ki,i-E K,i) are asymptotically normally 
distributed. The asymptotic variance of ^ X]r=i(^n,i ~ E V^„j) is given by 



Var i - E K.)) = - ^ Var(K,i) 

1 /i2/i2A 2(-l/WVul _Ll-A N 1/^1 



2 



n TT^yh^ 2 ij, In X n 

by Lemma O and . 

It is easy to check that in all three cases the approximation error ()3.1H1 is of smaller order 
than the asymptotic standard deviation in the theorems, provided A > 1/3. Hence this error is 
indeed negligible. □ 

4 Proofs of the lemmas 
4.1 Proof of Lemma 12.11 

Note that, for m = 0, 1, ■ ■ ■ , and any < e < 1, 



The exponent 3A is fairly arbitrary but it suffices for our purposes. By Lemma 3.5 of Van Es 
and Uh (2001) we have 

A 
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By applying integration by parts twice we get, for < e < 1, 



Hence 



A 



A 
A 



^m-X+l^{s/h)^/li 



'-(m-A + 1) I s'"-V^/'')'/'^rfs) 



^m~2\+l^(s/h)^/^l 



^ /,Ag(i//.)VM _ ^(^ _ A + l)/i2Ae(i/'^)V/^ 
A A"' 

+ ^{m - A + l)(m - 2A + l)h'^e^yh'^'' '^^ + o(/i3^e(i/'^)'/'^). 
A"* 



A A^ 



+ ^(^ _ A + l)(m - 2A + + o(/i3\(i/'^)'/^). 

A'^ 

This expansion is used repeatedly in the remainder of the proof. 
For m = we get 



cih) 



A A" 



A'^ 



which proves ()2.5|) . 

Furthermore, using (1 + x)^^ = 1 — x + + o(x^) for x | 0, we have 

J_ = + ^(A - l)h^ + ^(1 - A)(l - 2A)/i2^ + o(/ 

c(/;.) yU V A A^ 



A 
— / 



;^-Ag-(iM)VWl _ ^(A _ _ _ A)(i _ 2\)h 



A 



,2A 



A 



A2 



+ (^(A - 1)/.^ + ^(1 - A)(l - 2\)h'>^) + o(/.2^) 

A/,-Ae-(WV^ f 1 - ^(A - 1)/^^ + ^(A - A^)/^'' + o{h 
u V A A^ 
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Hence, by fj4.1|) for m = 1 and ()4.3|) . 



ES 



1 ' se^'/'^^'/^ds 



c{h) 



(l - ^(A - l)h' + ^(A - X')h'' + o{h'')) (l - ^(2 - X)h> 

+^(2-A)(2-2A)/i2^ + o(/i^)) 
l-^h' + ^{2-2\)e' + o{h''), 



which proves ()2.6|) . 

Similarly, by (j4.H) for m = 2 and (j4.3p . 

"1 



cih) Jo 

= (^i-i^{X-l)h' + ^^{X-X')h'' + oih' 

2 

(l - ^(3 - A)/i^ + ^(3 - A)(3 - 2X)h^^ + o(/i2^)) 

= l_2^/i^ + ^(6-4A)/i2^ + o(/i2^). 
A A^ 

Finally, by and (jO)) we get 

VaiS = ES^ -{ES)^ 
= 1 - 2^h' + ^(6 - AX)h'' + oih'') - (l - ^/i^ + ^(2 - 2A)/i^^ + o{h 



(4.4) 



(4.5) 



2A^ 



1 - ^h'' + 0{h 



2A\ 



which proves ()2.7p . 



□ 



4.2 Proof of Lemma 13.2 



Let e„ = — /i'^/^/ log/i denote a sequence of positive (for h < 1) numbers converging to zero. 
Note that for \t\ small enough we have 



A 



<A|A-l||t|. 



With t = fih'^v/X, for —e^Xh ^/ fi < v < 0, we have — < t < 0, and for n large enough 

A 1 // /, , nA 



((l + ^/i^t;) -l) -a| < A|A-l|^/i"|i; 
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and 



1 1 ^A 



<\X-l\-h-\l = o{l). 
/i 



< |A-l|^/iV 
A 



This implies that for these values of v, and n large enough, 



\fEAv)e-'' - 1\ 



X c{h) 
- A c(/i) 



1 1 



/i /i-^ 
1 1 



exp - — 1 + ^h^v -1) -v] -1 



A 



^'^ xcihy 

- A c(/i) 



1 1 



A c(/i) A 

< 3^/iV + 0(/i^), 
A 

where the remainder terms do not depend on v. 
Now note that 



0{h' 



and, by Lemma f2.H using (1 + s)^ — 1 < (A A l)s for —1 < s < 0, 



-e„{\/fj.)h 



(A^)h- 



H he 



^ c{h) J-{x/ti)h-^ 



,l/{,ih^){il+{t^/X)h^v)^~l)^^ 



- h-\l + 0(1)) / e'/^^'"^^^'+'^"-''>ds < - h~\l + o(l)) / e'/^^''"^^''''>ds 

J-1 J-1 



< 



X 

A A 1 



(l + o(l))e 



-e„(A/^){AAl)/i- 



0{h'). 
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Because the absolute value of the cosine is bounded by one we get 



E 



s„ cos h^-'E„{Xj - x)) - E^cos ( - ^ h^-'E{Xj - x)) 

< / \fEM-e''Ii-oo,oMdv 

/ \fEAv)-e''\dv+ / e^dv 

J -{\/lJ,)h->' J-oo 

j—en{X/ ix)h-^ 

fEAv)e-^ - l\e^dv + / IJeM ' e^\dv + e'^"/'')'^" 

J-{\/lx)h-^ 

< f (^h\^ + 0{h^))e''dv + 0{h^) = 0{h^), 



which proves the first statement of the lemma. 

The second statement can be proved similarly. □ 
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